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Abstract 



We suggest a method of constructing gauge invariant quark and gluon 
distributions that describe an abstract QCD observable and apply this 
method to analyze angular momentum of a hadron. In addition to the 
known quark and gluon polarized structure functions, we obtain gauge 
invariant distributions for quark and gluon orbital angular momenta, 
and consider some basic properties of these distributions and their mo- 
ments. 
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I. INTRODUCTION 



Since the initial realization that quark spin accounts for only a small fraction of 
the nucleon spin , theorists have struggled to define @ and suggest ways to measure 
the other components of the angular momentum in QCD. At first even the meaning 
of the quark spin contribution, E, to the nucleon spin was the subject of considerable 
debate. Now it is gradually agreed that different interpretations of £ reduce to 
renormalization scheme dependence made particularly irksome by the role of the axial 
anomaly in this channel. The gluon spin fraction, AG, seems to be well defined in 
the context of renormalization group improved parton models '! , where the polarized 
gluon structure function Ag(xBj, Q 2 ) has the meaning @@ of a scale dependent gluon 
spin distribution over the Bjorken variable xbjQ so that 

AG(Q 2 ) = J dx B] Ag{x Bh Q 2 ) . (1) 

Of course, Ag and AG are physical, gauge invariant quantities and experiments are 
being pursued to measure them. Nevertheless, there exists no local gauge invariant 
operator that could represent the space density of the gluon spin, a fact that has 
complicated the description of the gluon spin fraction. 

One may expect to incorporate the existing quark and gluon spin fractions into 
the nucleon angular momentum sum rule @ 

l - Y> + L q + AG + L g = ^ . (2) 

However, the situation with respect to the quark and gluon orbital angular momenta 
L q and L g is less satisfactory. There are no generally agreed upon definitions of L q 
and L g , and no associated gauge invariant definitions of their parton distributions @, 
/l 9 (^Bj,Q 2 ) and fL g (xBj,Q 2 )- Some authors use L q and L g given by the gauge 
variant decomposition of ref. ||, which appears in eq. (|5^ ) of our paper, in the gauge 
A + = 0, but alternative gauge fixing methods might also be considered. Besides, 
A + = does not completely fix the gauge, and L q or L g obtained this way are not 
separately invariant under the residual gauge transformation, that is not permissible 
for physical objects. 

We give a positive resolution to these questions in Sec. 3 of this paper, where we 
obtain scale dependent but otherwise unambiguous distribution functions for every 
term in eq. @ in a single framework. This will be contrasted with the conventional 
gauge invariant operator description of angular momentum in a gauge theory, 
where the decomposition (0) and the polarized gluon distribution Ag are obscure 



1 As usual, XBj = Q 2 /2P-q where P is the nucleon 4-momentum, q is the momentum 
transfer to the nucleon, and Q 2 = —q 2 . 
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and might even seem unphysical. But first we find it useful to develop a general 
approach to the construction of parton distributions associated with an arbitrary 
physical observable that may characterize quark and gluon fields. This approach is 
presented in the Sec. 2. The premises and results of Sees. 2 and 3 are summarized in 
the conclusion. 

II. STRUCTURE FUNCTIONS FOR A PHYSICAL OBSERVABLE 

We want to describe some additive physical quantity T relative to a deep inelastic 
process, such as inclusive leptoproduction, where the Bjorken variable XBj = Q 2 /2P-q 
is fixed by the kinematics of the experiment (scattering angles and energies) and the 
information about the hadron state is limited to forward hadron matrix elements of 
certain operators. In a naive free parton model and the Bjorken limit q~ — > oo 

ZBj = p+, (3) 

where p + and P + are the "+" components of the struck parton and the target mo- 
menta correspondingly. Therefore a similar experiment sensitive to T would directly 
provide "the density of T" carried by quark and gluon partons with the momentum 
p + = XBj-P + , rather than T itself. With this in mind, we construct gauge and boost 
invariant expressions for such densities that we call quark and gluon T distributions, 
fr«( x Bj) and /rs The observable T can be defined by any hermitian operator 

constructed from quark and gluon fields or as a generator of some fields transforma- 
tion. We do not associate T with a specific experiment measuring it. Therefore, in 
this paper we will not give a universal answer to the question: In what experiment 
one should measure a general observable, but we show what knowledge about T can 
in principle be obtained in a deep inelastic process. 

In order to construct the parton distribution of F in QCD we should make clear 
what a quark or gluon parton is in an interacting theory. The subtlety here is in 
the issue of gauge invariance: a pure quark field in one gauge is a superposition of 
quarks and gluons in another. Different ways of gluon field gauge fixing predeter- 
mine different decomposition of the coupled quark-gluon fields into quark and gluon 
degrees of freedom. Similarly, one can generalize a gauge variant nonlocal operator, 
obtained, for example, by Fock space decompositions of free fields, to more than one 
gauge invariant expressions, raising the problem of deciding which is the "true" one. 
We choose to define quark and gluon partons in fully interacting QCD such that their 
distributions in XBj are given by the unambiguous experimentally accessible struc- 
ture functions '§ q(xBj) and g(xBj)- We find that this definition also reproduces the 
measurable polarized structure functions Ag(a;Bj) and Ag(xBj)- 
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The required definition of QCD partons is also very appealing from mathematical 
point of view. The "gauge covariant" quanta of quark and gluon fields with a def- 
inite momentum p + = xr- s P + are the eigenstates of the generators of the covariant 
translations in the x~ direction: 

Tl(a~) : ip(x) — > U(x, x + a~) ip(x + a") , , , 

Ti(ar) : D x {x) -> U(x,x + aT) D x (x + a~) U(x + a~,x) , ^ ' 

where x + a~ is the shorthand for (x + ,x~+ a - ,^ 1 ,^ 2 ), U(x, x+a~) and U(x+a~, x) 
are the Wilson lines along the straight path between x and x + a", 

U(x, x + a~) = Pexp[ig [* d£-A + {£)] (5) 

with a similar equation for U(x+a~ ,x) = W(x,x+a~), and D\(x) = d\ — igA\(x), 
so that the second line in eq. (||) defines a transformation of the gluon field A\(x). 
The transformation Tl commutes with Tl and each of them forms a one-parameter 
Lie group: 

[T?(a-),2?(6-)] = 0, 
Tl{a^)Tl{al)=Tl{al + a 2 ) , (6) 
T^b 2 )T»_(b- 1 )=Tl(b- 1 +b 2 ) . 

These properties are easily established from the fact that T£ (a~) leaves A + (x), and 
consequently the Wilson line (|5]), invariant, as evident from the infinitesimal form of 
the transformations (fD 

Slip = a'D-ip , , s 

S 9 _A x = a-F_ x S g _A + = 0. [) 

Obviously the transformations Tl and Tf, as defined by eq. (Q) or eq. (0), commute 
with gauge transformations. 

Action of both Tl(a~) on the quark field and T£ (a~ ) , with the same parameter a~ , 
on the gluon field yields ordinary translation of these fields by a - plus an overall 
gauge transformation. This additional gauge transformation does not affect any gauge 
invariant observable, therefore the gauge covariant translation of every field in the 
system is physically equivalent to their usual translation. This is not true, however, 
when one considers the transformations (|j) for quarks and gluons individually. For 
example, application of Tl(a~) to the quark field ip(x) without changing the gauge 
field A\{x) gives a physically different field configuration from what would be obtained 
by the usual translation if)(x) — > ip(x+a~), because a change of the i/j(x) phase without 
a compensatory change in A\(x) is not a symmetry of QCD. 

The quark and gluon fields if)(x) and A\{x) can be decomposed into a sum of 
one-dimensional irreducible representations of the groups Tl and T£, which are the 
eigenstates of the corresponding generators. We will see that this decomposition 
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specifies the desired parton decomposition in an interacting theory. The irreducible 
representations of Tl and T£ become trivial in the light-cone gauge A + = where 
these transformations coincide with ordinary translations. Let us denote fields in this 
gauge by the superscript l LC\ Then 

and the decomposition that we are seeking for is simply 

^ LC (x) = y _ e -** + »-^(fc+ i) , (9) 



where 



Ax {X)= I 17 6 Mk + ,x) (10) 



x = (x , X , X ) . (11) 

Since Tl and TO. commute with gauge transformation, their action on fields in an 
arbitrary gauge 

Mx) = e ia{x) i) LC (x) , A x (x) = e ia{x) (A L x c (x) + - d x )e- ia(x) , (12) 

9 

amounts to translation of the "physical" components ip LC and A x c according to 
eq. (|]) with the unphysical phase a(x) left invariant: 

Tl(a~),Tl(a~) : a(x) — > a(x) . (13) 

Therefore, the "parton decomposition" of fields in an arbitrary gauge is simply ob- 
tained by a gauge transformation of eqs. @ and (|l(J). These equations can be viewed 
as transition from x to (k + , x) representation which is more suitable for deep inelastic 
scattering, but this transition is done in a way respecting gauge invariance and coin- 
cides with brute-force Fourier transformation along x~ direction only in the specific 
gauge A + = 0. 

Our use of the adjective "physical" for fields in the gauge A + = has been some- 
what premature. There remain gauge transformations, for which the phase a(x) in 
eq. flP2"| ) does not depend on x~, that preserve the condition A + = A_= 0. All true 
physical observables must be invariant under this residual gauge symmetry as well. 
The fields ip(k + ,x) and A\(k + ,x) defined by eqs. © and ( |10|) transform under it as 

ip(k + ,x) -> e ia(x ^(k+,x) , 

A x (k + ,x) -> e ia ^(A x (k + ,S:) + d x )e~ ia ^ . ^ 
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Notice that the inhomogeneous term in the residual gauge transformation of A\ ap- 
pears only at k + = 0. We find it useful to split the field A\(k + , x) into a part G\(k + , x), 
transforming homogeneously, and a gauge field of the residual gauge group, A\(x): 

A x (k + , x) = 2n5(k + )A x (i) + G x (k + , x) , (15) 

so that under a residual gauge transformation 

A x (x) -> e ia ^(Ax(x) + - dx)e~ ia ^ , G x (k + ,x) -> e ia ^G x (k + , x)e~ ia ^ . (16) 

The decomposition ( |T5j ) becomes unambiguous if we require 

Gx(k + ,x)\ k+=Q = 0. (17) 
Then combining eqs. fllCf ) and (0) one obtains 

^a(x) = / dz~ Af?{z) . (18) 

The fields ip(k + , x), Gx(k + , x), and the residual gauge field A\(x) provide the basis in 
the field space that we use to construct gauge invariant quark and gluon distributions 
for physical observables. 

We associate a physical quantity T with the field transformation generated by the 
quantum operator corresponding to T. We do not insist on the invariance of the QCD 
action under this field transformation. If it is a symmetry of the theory then there 
exists a corresponding conserved Noether current, and its charge, which will be the 
operator of T, commutes with the Hamiltonian of the system. If it is not, the quantum 
expectation value of T for a particular state is still a well defined quantity, and spin 
or orbital angular momentum, which are not separately conserved in a relativistic 
system, provide such examples. 

Of the possible choices of T, only those that commute with the parton-defining 
transformations Tl and T£ can be interpreted as physical properties of quarks and 
gluons distributed in x B j ~ that is, as "parton distributions of T" . Thus we restrict 
ourselves to T's that commute with Tl and : 

[6'L,6 r ] = [6 9 _,6 r ] = 0. (19) 

This allows complete description of F in a hadron in terms of its distributions over 
partons, /r^s^Bj)- 

These distributions, i.e. parton densities of T, are easily constructed using the 
following trick: We apply T independently to each subspace of fields with a defi- 
nite k + . Remembering that by eq. (|19"D T can be diagonalized in the basis of Tl and 
Tl eigenf unctions, we consider the set of independent transformations labeled by a 
continuous parameter XBf 
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k+- 



-p 



p+ 
k+- 



x Bi P+ 



(20) 



I'Bj 



P+ 



and 



r(^Bj) a 

Oris -A\ 







(21) 



where T q and V 9 are the generators of 5r restricted to the fields if)(k + , x) and G\(k + , x) 
with k + = x-Q i P + . The generators T q and V 9 can be operators acting on spinor, 
Lorentz, color, etc. indices as well as on the coordinate x _L =(x 1 ,x 2 ). For simplicity 
we slightly limit our choice of V further, assuming that T q and V 9 do not depend on 
the parameter q + and do not involve the coordinate x + or derivatives with respect to 
it. This is not absolutely necessary but it simplifies the following equations. We will 
not need such complications for any of the quantities that are required for angular 
momentum description. 

The generators of the field transformations defined by eqs. (^) and fl2~l|) are asso- 
ciated with a physical observable which is their normalized expectation in a hadron 
state \P). These observables will be shown to be the V parton distributions that we 
are looking for. Applying Fourier transform to eqs. (EDf) and Q2"T|), one finds the action 



of (5p XBj ' ) on the conventional fields ip LC (x) and A x c ( 



x 



5 { * q B ^ LC \x) 



ieP^ 



ieP^ 
~27 



J d£-e iXBiP+ t T q ^ LC {x + r) , 

J dre iXBjP+r rf (a l x c (x + r ) - a 



These transformations are generated by the light-front charges 



Q« = J d 2 x ± J dx-J + (x) 



of the corresponding canonical currents^ : 

j&texBj) = C- / dre" BjP+r ^h / T^(x + r) 

ZTT J 



Jp 9 (x; x-q^ 



iP + 
~2T 



dCe lXBiP+r F^ x (x)T 9 x x \A x (x + T) - A x (x 



(22) 
(23) 

(24) 

(25) 
(26) 



2 By the canonical current for a field transformation 5r<& = ieT(&) we mean 

fields 
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It is necessary to use the light-front charges instead of the usual Q = J d 3 x J°(x) 
because the transformations fl2"2|)-(|2"3]) are non-local in x~ . We can imagine boost- 
ing ourselves to the infinite momentum frame, K', where the target hadron has an 
infinitely large momentum in the x' 3 direction: 

p> = (Ep, 0, 0, P' 3 ) with P' 3 -> +00 . (27) 

In this limit the integration plane in 

Q' = ( dVj'V) (28) 



"approaches" the On the other hand, at least for fields confined in a massive 

hadron, 

lim Q 1 = Q (+) , (29) 

therefore is the correct generator. One may also check this using the light-cone 
formulation of a gauge theory E§. The light-front charge is invariant under 
boosts along £3 and can be evaluated in any frame in which the hadron has zero 
transverse momentum, P 1 =P 2 = 0, including the rest frame. Notice that if the cur- 
rents (|25|), (|26| ) were conserved, their usual charges, Q, would be frame independent 
and equal to the light-front ones, Q^ + \ 

Now we define the quark or gluon distribution functions for the physical observable 
T in a covariantly normalized hadron state |P), (P'\P) = 2Ep(2n) 3 5 ( - 3 \P' — P), as 
the normalized expectation 

/r-M = (P| Q { r + Mx Bj ) |P) = 



1 



2P+(Jd 2 x ± Jdx 



(P| J d 2 x ± J dx-J+, g (x;x Bi ) |P) . (30) 



The distributions are normalized so that /r^^Bj) = Tp5(l— xbj) in a free parton 
model with |P) replaced by a quark or gluon parton which is an eigenstate of the 
transformation (|20| ) with the eigenvalue Tp. With such normalization, the sum of the 
first moments of the quark and gluon T distributions gives exactly the expectation 
value of the observable T for the hadron state |P) in the infinite momentum frame: 

1 A (f f \ ( P l^r + V) ,. (P'\Qr\P') ^ 

.! ^ (/r9 + /rs) = (P|P) = Ao (P>\P>) ■ (31) 



The objects defined by eq. ( |30[ ) are Lorentz boost invariant, provided T g and T 9 
transform as implied by their Lorentz and spinor indexes. They are invariant under 
residual gauge transformation if T q and V 9 transform homogeneously as: 
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(32) 



Some examples are T q = 1,7 M ,<9_|_ — igA± and = 5*; we will have more ex- 
amples in Sec. Ill . When /r«(£Bj) an d /rs(^Bj) are invariant under the residual 
gauge group, they are by construction the same in an arbitrary gauge. An explic- 
itly gauge covariant expression for the distributions is straightforwardly obtained by 
inserting the Wilson lines U(x, x+(,~) into eqs. (|25| ) and (|26| ) as ip. . .ip — > ipU. . .ip, 
F. . .A — > Ty[FU. . .AU']. Even though the currents in eq. fl3U| ) are not necessar- 
ily conserved and the charge operator depends on x + , the diagonal matrix element 
between 4- momentum eigenstates \P) are x + independent. 

The composite operators in eqs. (p5| ) and ([26|) have ultraviolet divergences @, 
which arise in perturbation theory as infinite loop integrals over k± in the products 
of two fields at the same x x . These divergencies are regulated by a cutoff at k\ ~ 
Q 2 = — q 2 leading to factors like a(/i 2 ) log (Q 2 / fi 2 ). In practice they are handled with 
standard methods of operator renormalization. Renormalization of objects similar to 
eq. ( p0|) is described in the work II by Collins and Soper. 

For our purposes we take the state \P) in the definition of distribution functions 
to be the hadron helicity eigenstate with the maximal possible helicity, e.g. +| for a 
nucleon. Of course, for spin independent structure functions one could equally well 
average over helicities. For some problems, for example transversity distributions @, 
one must consider off-diagonal matrix elements between states with different helicities 
but the same 4-momentum P. 

Notice that the very last term, A\(x), in eq. (|26|) does not contribute to the T 
distributions. Indeed, it produces the following term in /rs(^Bj) : 

iJQEBj) (P|/dV/dx^dy-F+ A (x)rfA x (y-,x) \P) 
f^)\A ~ 2P+ Jd 2 *M/d*-) 2 ' ( } 

where eq. ( pf) was used. Because of the xbj delta-function, this term trivially vanishes 
for XBj 7^ 0. One might worry that its value at Xbj = could affect the first moment 
of /re^Bj)- In fact, it does not: In the A + = gauge F +X (x) = d^A x (x) and A x (x) 
remains finite when x~—>±oo, thus the expression (|33|) vanishes as an infinite volume 
average of the derivative of a bounded function^. 

The familiar structure functions q(xBj), Ag(xBj), g(x B j), Ag(x B j) are instantly 
reproduced from eq. fl3U|). Indeed, the quark number and chirality are counted by 
T q = 1 and T q = 7 5 correspondingly. Translating the current matrix elements to 
the origin and canceling out the volume factors in the numerator and denominator of 
eq. (EH) we find 



3 In this discussion we treat gluon fields as if they where classical. Since eq. (33) deals 
with physical expectation values, the conclusion A/rs (xBj)|_4 = should remain valid for 
the quantum case. 
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H = 1 =► /.(xbj) = — ^ / dre^ P+r (P|4(0)^ + (D l-P) , (34) 

P = 7 5 =► /A 9 (x Bj ) = ^= / dT e^ P+r (P|4(0)7 5 ^ + (D l-P) (35) 

that are recognizable 111 as unpolarized, g a (xBj), and polarized, Ag a (a;Bj), quark dis- 
tributionsQ. One can easily write down the corresponding distributions q a {xBj) and 
Aq a (xBj) f° r a particular flavor a = u,d,... complementing r 9 's in eq. (|35|) with a 
flavor projector Pf), = 5f a 8p a . 

Similarly we proceed with the gluon structure functions g(xBj) and A^(xbj)- For 
the first one = 5* so that 

/,M = ^ / dr e^ P+r (P\F + \0)A x (C) \P) • (36) 

Substituting e JXBjP+ ^ by (ixBj-P + ) _1 g|^e J:!;BjP+ ^ , integrating by parts, and remem- 
bering that d-A{ c = F\ LC = -F X +LC we obtain 

= / dCe^ p+ Z-(P\F+\0)F x +(C) l-P) + [(• • Ol£=S] • (37) 

The first term on the right hand side of eq. ( |3"7| ) is a standard form of the unpolarized 
gluon distribution It has been argued H that the second, surface, term in 

eq. ( |37|) can be dropped. We postpone the discussion of the gluon spin dependent 
distribution Ag(xBj) until the next section, where it appears in a natural way. 



III. ANGULAR MOMENTUM IN DEEP INELASTIC SCATTERING 

The application of our formalism for obtaining ^Bj-distributions that describe an- 
gular momentum in a hadron is straightforward: one should identify the correspond- 
ing transformations of the quark and gluon fields tp and G and substitute their gener- 
ators T 9 and V 9 into the master formulas (|3"0|). As mentioned earlier, we are interested 
in transformations that are diagonalized by the "parton decomposition" (0), ( |10|) or 
equivalently that commute with the gauge covariant translations given by eq. (f|). 
Among three-dimensional rotations these are the rotations about the x 3 -axis gener- 
ated by J 12 that transform the fields ipip+^x) and G\{p + ,x) as 

W = ~e [i- 1 ^ 12 } + {x 1 d 2 - x^}] $ , 
5jG x = -e[{-(5l5X-5l5X)} + {x 1 d 2 -x 2 d 1 }]G x . 1 ] 



4 In eq. (|35| ) and later = ^7 j + ip is the "good" component of the Dirac field ip, which 
is the dynamical field in light-cone quantization @ . 
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We consider the transformations generated by each of the four terms surrounded by 
braces in eq. (|38|) independently of the others. The first generators in each of these 
equations are naturally associated with the field spins whereas the second are the 
standard generators of orbital angular momentum. In order obtain to gauge invariant 
structure functions we must maintain covariance with respect to the residual gauge 
transformations, in which the local gauge transformation parameter does not depend 
on x~ . This can be achieved replacing <9j by the residual gauge covariant derivative 

T>i = di- igAi . (39) 

Similar to the situation at the beginning of the previous section, the replacement (|39| ) 
adds only an overall unphysical gauge transformation to the rotation ( |38| ) when 5j 
applies to both quark and gluon fields with the same infinitesimal angle e. Also 
notice that by construction of A\ in eq. (|T8| ) we have T>i = di in any gauge where 
lim x -_ ±00 Af c (x) = 0. 

Splitting the transformation ( |3lf ) into the four independent ones as described, we 
define 

E : <5 s t/> 

L q : 5l$ 
AG : 5agG\ 
L g : h g G x 

The corresponding generators are 

T s = a 12 = 7 °7 3 7 5 , r Lq = (xHV 2 - xHV x ) , 

(r AG )A x = -ie+~x x , (r Lg )x x = (xHv 2 - x 2 w x )8l , 

where £^ ua f 3 is the antisymmetric Levi-Civita tensor with £ 0123 = +1 and e + ~ 12 = —1 . 

Following tradition the generator for E is defined without the \ factor that ap- 
pears in front of a 12 in eq. fl3~8l) . Substituting these T's into eqs. (f25|)- (|26|) , (|30|) and 
remembering that the A\(x) term in eq. ( p6|) can be dropped, we obtain the following 
structure functions for the quark and gluon spins and angular momenta:^ 



tea 



12 



ip = ie'y 7 7 ip 
-e{x x V 2 - x 2 Vi)t]j ; 
e{%% - 5 2 5f)G x = 
-t[x x V 2 - x 2 V x )G 



e 6 + - x xG 



(40) 



5 We remind the reader that = ^7~7 +, (/', the hadron \P) is polarized along its momen- 
tum P, and the fields tp, A\, and F^ u in the structure functions below are either connected 
by a Wilson line to make these expressions gauge invariant or are evaluated in the gauge 
A + = 0, where the Wilson line is unity. 
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/s(XBj) = 2^2 / d ^ e " BjP+r <Wl(0)7 5 <Mr) l-P) ; (42) 

h ^ ) = 2WHf¥^) ; (43) 

/ag(x Bj ) = dre l ^ P+r (P|F +A (0) £ + - A M x (r) l-P) ; (44) 
, , , z/dC~e^ p+ ^(P|/dVF +A (rr ± )(x 1 zX> 2 -x 2 zPi)A A (a; ± + C~) l-P) ,„ K x 

/Ls(XBj) " 47T (JdV) • 

In these expressions we have translated fields to the origin and canceled overall inte- 
grals where it is possible. One can replace the vector potential A\ in eqs. (|44]) and (|45| ) 
in favor of the gluon field strength F x + as was done when we derived eq. (|37|). Drop- 
ping the surface terms, we obtain: (below F*^ v = ^e^ ua ^F a p) 

/ag(x Bj ) = / dre^ P+ ^(P|P +A (0)P*/(r) IP) , (46) 

fd£-e ix vi p+ S-(P\Jd 2 x ± F +x (x ± )(x 1 iV 2 -x 2 iV 1 )F x + (x ± +£-)\P} 
/L * (XBj) " 47TXBjP+(/dV) ' (47) 

These are our principal results. Equations (^)-(|47|) give a complete set of distri- 
bution functions needed to describe a hadron's angular momentum. These expressions 
are automatically gauge and boost invariant as a consequence of our method. Our 
expressions for orbital angular momentum differ from the naive gauge invariant form 
in which Di appears instead of ZV This distinction is essential because T>i depends 
only on x whereas Di is also a function of x~, which is impossible to interpret in a 
parton model formalism. Our distribution functions can be written in an explicitly 
gauge-invariant form with use of the Wilson connection. We devote the rest of the 
section to analyzing these expressions and comparing them with earlier works. 

In a non-interacting parton model the structure functions /s, fi, q , /ag, and fi g 
describe the ^Bj-distributions of spin and orbital angular momentum x 3 -projections 
carried by quark and gluon partons respectively. This is also true in the interacting 
theory when one defines quark and gluon partons as excitations of the fields ip{p + , x) 
and G\(p + ,x) and their spin and orbital momentum as the charges of the transfor- 
mations (f40|). We demonstrated that these are the definitions that reproduce the 
familiar unpolarized and polarized quark and gluon structure functions. According 
to eq. (PH), the first moments of the structure functions 

l 

i dx B j /r(^Bj) , (48) 
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should give the total quark/gluon spin/orbital momentum in the hadron in the infinite 
momentum frame. The sum of the first moments obeys the angular momentum sum 
rule in eq. (|2|), where 1/2 on the right hand side is the ^-projection of the nucleon 
spin and should be replaced correspondingly for a hadron with a different spin. The 
first moment of /g 

S = ^ T (P|^l(0)7 5 ^ + (0)|P) , (49) 

indeed equals the total x 3 -projection of the quark spin 

(P'\ J d 3 x 4){x)^i)(x)\P') 1 
(P'\P') ~ 2E P 



as P' 3 — > +oo. 

The situation for L q , AG, and L g is a little more subtle. The first moments of the 
corresponding structure functions are 

L q = ^ 1 d 2 x±) ( p \ J d 2 x ± ^ + (x ± )(xHV 2 - xHV^x^P) , (51) 
AG = ^(P\F +x (0)e + - x xA x (0) \P) , (52) 

L g = 2p+ » (P[ J d^F+\x^)(xHV 2 - xHV^A^) \P) ■ (53) 

Up to the replacement di <-> T>i, the eqs. fl49]) and d5l|) -(|53|) reproduce the normalized 
forward matrix elements of the four operators in the gauge variant angular momentum 
decomposition by Jaffe and Manohar % 

J 3 = J d 3 x { -Vi7 3 7 5 -0 + V f (x x idftp + {Ex A) 3 - E k (x x d) 3 A k ] (54) 

evaluated in the A + = gauge and the infinite momentum frame. This result is quite 
welcome in view of previous studies of angular momentum in deep inelastic 

scattering. This implies that many previous calculations, including Q 2 evolution and 
asymptotic behavior §h@, obtained in A + = gauge are transferable to our case with 
only minor modifications to incorporate ck — > ZV 

On the other hand, it is also known til that the total angular momentum^ 



6 In the following eq. ((5^) Q^ u stands for the symmetric energy-momentum Belifante 
tensor. 
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can be decomposed into three gauge- invariant vectors: 

J = J d 3 x I ^77 5 ^ + i)\x x iD)ip + x x (E x B)) , (56) 

where 

E k = F k0 ^ B k = _}_ e klmp im ^ 

are the color-electric and color-magnetic field vectors. The second term in eq. ( EBp is 
an alternative candidate for a definition of the quark orbital angular momentum. The 
third, gluon, term does not allow any further separation into local gauge-invariant 
operators that would describe spin and orbital momentum of gluons^. The d 3 x inte- 
grals of each of the three terms in eq. (|56|) are still not covariant under Lorentz boosts 
so one have to specify a frame in which they are evaluated. These quantities are in 
principle measurable and can be related to form-factors of deeply virtual Compton 
scattering @ . Nevertheless, the decomposition fl5*rj|) is not favorable for deep inelastic 
scattering experimental conditions in the sense that the separate terms in eq. (|56|) mix 
partons with different p + and thus can not be presented in terms of parton densities. 
Of course the greatest shortcoming of the decomposition of eq. (^) is its inability 
to describe the gluon spin spin distribution Ag(xBj) which is clearly gauge invariant 
and can be measured in a variety of processes. There is no such a difficulty in the 
suggested approach where Ag(xBj) = /ag(^Bj)- 



IV. CONCLUSION 

Typically in (semi)-inclusive deep inelastic scattering the Bjorken variable xbj is 
fixed, and only forward hadron matrix elements are accessible. In order to accommo- 
date these experimental conditions, we propose to describe a hadron in terms of those 
observables V that are diagonal in the basis formed by quark and gluon partons. We 
then construct gauge invariant xbj -distributions associated with such an observable. 
The gauge covariant definition of parton states in fully interacting QCD are chosen 
so that the parton distributions are given by the physical structure functions q{xBj) 
and ^(^Bj). To satisfy this requirement, quark and gluon partons must be eigenstates 
of the generators of covariant translations Tl and T£ , eq. (^), and the T distributions 
are given by the formulas 



We were able to give gauge invariant definitions for the gluon spin and angular momentum 
projections only on the axis x 3 defined by the external momentum q. It is not possible 
to define gauge invariant tensors of gluon spin and orbital momentum independent of a 
specified direction in space-time. 
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(-P|/dV-/dar/d£- e ix ^ p+ ^ ^{x^+T^x + f") |P) 
4tt(/ dV / dar 



/r«(^Bj) = 7Z7732-I ( 5 * 



and 



(P|/d 2 x^/dx"/dre-B^ + r F+ A (x)rfFx+(x + r) )p) 
/r9(XBj) = 47rx Bj P+(/dV/dx-) (59) 

in A + = gauge. They may be put in a gauge- invariant form by trivial use of Wilson 
lines. In these formulas the generators T q and T 9 transform under the residual gauge 
group as given by eq. fl52p and are uniquely related to the original quantity T. 

A hadron angular momentum is completely described by the four scale dependent 
XBj-distributions given by eqs. (P^)-(]4"3|), (|5)-(H7[). Two of them, / E and /ag> coincide 



with the polarized quark and gluon structure functions Aq and Ag. The other two 
are naturally regarded as the xbj distributions of quark and gluon orbital angular 
momentum. They are well defined physical objects and are gauge invariant by means 



of the residual gauge covariant derivative P, = di — igAi, where A is given by eq. fllSP . 
They are observable in principle and can be calculated in models or within lattice 
QCD. However, it remains an open question whether the XBj-distributions of quark 
and gluon orbital angular momentum which we have defined can be measured in a 
practical experimental process. 
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